% 计算辐射函数 F
function F_theta = F(Theta)
    b = 2;  % Boresight gain参数
    if Theta >= 0 && Theta <= pi/2
        F_theta = 2 * (b + 1) * cos(Theta)^b;
    else
        F_theta = 0;
    end
end

